Complex shape evolution of electromigration-driven single-layer islands 
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The shape evolution of two-dimensional islands through periphery diffusion biased by an electro- 
migration force is studied numerically using a continuum approach. We show that the introduction 
of crystal anisotropy in the mobility of edge atoms induces a rich variety of migration modes, which 
include oscillatory and irregular behavior. A phase diagram in the plane of anisotropy strength and 
island size is constructed. The oscillatory motion can be understood in terms of stable facets which 
develop on one side of the island and which the island then slides past. The facet orientations are 
determined analytically. 
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The manipulation of nanostructures by macroscopic 
forces is likely to become a key ingredient in many nan- 
otechnology applications. Understanding the influence of 
external fields on the shape evolution of nanoscale sur- 
face features is therefore of considerable importance. As 
a first step in this direction we analyze here the effects of 
an electric current on single-layer islands on a crystalline 
surface. The islands evolve under surface electromigra- 
tion, the directed motion of adsorbed atoms due to the 
slight force transmitted by collisions with the conduction 
electrons in the sample 1] . 

Electromigration along interfaces and grain boundaries 
is the most persistent and menacing reliability problem 
in integrated circuit technology 0, El ■ Correspondingly, 
much work has been devoted to electromigration-induced 
void formation and breakdown in metallic conductor lines 
0, and the capacity for quantitative numerical model- 
ing has been demonstrated at least for simple void ge- 
ometries 0, 0- A major obstacle to achieving predic- 
tive power in such studies, however, is the insufficient 
control over the complex internal structure of the poly- 
crystalline samples. Hence an important motivation for 
investigating electromigration-induced effects on simple, 
well-controlled nanoscale morphologies, such as stepjpat- 
terns on vicinal surfaces and single layer islands , is 
to bridge the gap between the microscopic mechanisms 
of electromigration and their consequences on technolog- 
ically relevant length and time scales. 

Electromigration of islands has been modeled previ- 
ously using Monte Carlo simulations Q and continuum 
theory |9( . The continuum approach to island shape evo- 
lution, which treats the island edge as a smooth curve, 
has been successfully applied to a range of problems in- 
cluding the diffusion jlOj and sintering 



1 1 1 1| of islands, and 
the pinch-off of vacancy clusters 12]. Here we focus on 
the regime of periphery diffusion (PD), where the domi- 
nant kinetic process is the migration of atoms along the 
island boundary. The shape then follows a local evolu- 
tion law, without coupling to the adatom concentration 



on the surrounding terrace. 

We extend the model of by including crystal 
anisotropy in the adatom mobility. It was observed re- 
cently in the context of step flow growth [13( that crys- 
talline anisotropy can change the behavior of step pat- 
terns in a qualitative way. In the present case, it leads to 
the unfolding of a remarkable richness of dynamic phe- 
nomena: In addition to the scenarios of steady motion 
and island breakup [3, |9| observed in previous work, we 
find spontaneous symmetry breaking, oscillatory shape 
evolution, and complex migration trajectories where dif- 
ferent modes of motion alternate in a periodic or irreg- 
ular fashion. This highlights the importance of properly 
including anisotropy in the modeling of boundary evolu- 
tions. Oscillatory shape dynamics has been seen previ- 
ously in a numerical study of void electromigration |l4| , 
and transitions to quasiperiodicity and chaos are known 
to occur in directional solidification LLij ■ To the best 
of our knowledge, however, our work provides the first 
example of complex shape evolution for a closed contour 
subject to purely local dynamics. 

In the PD regime, the normal velocity v n of the island 
boundary satisfies the continuity equation 
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Here s denotes the arclength along the island edge, and 
O the atomic area. The square bracket multiplied by the 
edge atom mobility a is the total mass current along the 
boundary, which is driven by the tangential derivative 
of the chemical potential A/i = fl'jK and the electromi- 
gration force q*E t ; 7 is the edge stiffness, k the local 
curvature, q* the effective charge of an edge atom, and 
E t the tangential component of the local electric field. 
The crystal anisotropy of the surface enters through the 
dependence of 7 [13 an d & HH El on the orientation 
angle 8 of the island edge. 

For atomic layer height islands on the surface of a thick 
sample, the island boundary has a negligible effect on the 
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FIG. 1: Phase diagram of migration modes in the plane 
of anisotropy strength S and island radius Ro, for sixfold 
anisotropy (n = 6) and the field aligned with a direction of 
maximal mobility (a = 0). For each point on a grid of resolu- 
tion 0.5 x 0.5, the evolution of the island was followed until the 
asymptotic mode could be identified. We distinguish between 
straight stationary (ss), oblique stationary (os), oblique oscil- 
latory (oo), zigzag (zz) and complex oscillatory (co) motion. 
In the region bu islands break up. The cross on the i?o-axis 
indicates the linear instability of the circular solution in the 
isotropic case. 



electric field; this is in contrast to the modeling of voids 
in conductor lines, where the coupling of the void shape 
to the electric field leads to a manifestly nonlocal prob- 
lem 0,0, El- Here we can take the field to be of constant 
strength Eq and aligned along the x-axis. Letting 9 de- 
note the angle between the normal of the island edge and 
the y-axis (counted positive in the clockwise direction), 
this implies E t = Eq cos(0). 

Together with the specification of 7 and a, to be ad- 
dressed below, this completes the definition of the local 
boundary evolution Q). Comparing the two terms inside 
the square brackets, we extract the characteristic length 



scale [iilSlliSl Ie = /\q*Eo\, which gauges the rela- 
tive importance of capillary and electromigration forces; 
electromigration dominates on scales large compared to 
Ie- Below all lengths are reported in units of Ie- 

The isotropic version of Q, with 7,17 = const., 
has been studied previously by Suo and collaborators 
EiJ ESI El- A circular island moving at constant ve- 
locity is stable for (dimcnsionless) radii R < R c « 3.26 
|2~lI | . Beyond the instability a bifurcation to two branches 
of non-circular stationary solutions occurs |2fJ . Numeri- 
cal integration of the time-dependent problem shows 
that only one of the branches, corresponding to islands 
elongated in the field direction, is realized. At large radii 
island breakup occurs, mediated by the outgrowth of a 
finger of the kind found in 0] . 

We now turn to the effects of crystal anisotropy. 
Throughout this paper only the anisotropy of the adatom 
mobility a will be taken into account, while the edge stiff- 
ness 7 is kept isotropic. This is motivated partly by the 
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FIG. 2: Stationary shapes for S = 2 near the transition from 
ss to os behavior. Arrows indicate the direction of motion. 



fact that the anisotropy in a is found experimentally (to 
the extent that it has been investigated) to much exceed 
that of 7 [lfij . and partly by our desire to disentangle 
kinetic (a) and thermodynamic (7) effects j^. For the 
kinetic anisotropy we employ the functional form [j| 

<t{9) = a max (l + Sy^l + Scos 2 [n(6 + a)/2]}. (2) 

Since the prefactor er max only sets the time scale, the 
relevant parameters in @ are the anisotropy strength S, 
the number of symmetry axes n, and the angle a between 
the symmetry axes and the electric field direction. The 
natural time unit is ts = 1% / '(o"max7^ 2 ) 3- 

The simplest solutions in the anisotropic case are sta- 
tionary islands moving at constant speed, which satisfy 
the equation v n = V sin (9 + 4>); here the angle (f> accounts 
for the fact that the island does not necessarily move in 
the field direction. A complete analysis of stationary is- 
land shapes has been achieved in the limiting case of zero 
stiffness j2^. For an even number n of symmetry axes 
smooth stationary shapes are found for small S, while for 
larger anisotropy the shapes develop self-intersections; no 
stationary shapes exist when n is odd. The migration di- 
rection generally lies between the field direction and the 
symmetry axis of the anisotropy. 

Despite their mathematical interest, these results are 
of limited applicability to real islands, because all sta- 
tionary shapes are wildly unstable when 7 = 0. In the 
remainder of the article we therefore focus on the numer- 
ical solution of the full, time-dependent problem with 
7 > and a(9) given by @. Two complementary nu- 
merical algorithms have been employed. For relatively 
small islands a finite difference scheme described in |5j 
was found to be most efficient, while for large islands we 
rely on the better stability properties of a semi-implicit 
adaptive finite element algorithm. [24^ . The full mutual 
consistency of the two approaches has been checked. 

Most results have been obtained for n — 6 and a = 0. 
This leaves the anisotropy strength S and the initial con- 
dition for the deterministic shape evolution to be spec- 
ified. Extensive calculations show that the dependence 
on the precise initial shape is minor |25| . and hence the 
initial condition can be characterized by the radius i?o of 
a circular island of the same area; in practice, we usually 
start the calculation from a slightly distorted circle. 

The phase diagram in FigQ gives an overview of the 
observed migration modes in the S — i?o-plane ■ For 
small islands (Rq < 2) the evolution converges to a sta- 
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FIG. 3: Snapshots of oo motion for Rq = 4 and S = 1 (upper 
panel) and zz motion for Rq = 3.5, 5* = 0.5 (lower panel), 
taken at time intervals At = 20. In both cases the perimeter 
displays simple oscillations, as in the bottom panel of Fig|l] 

tionary shape which moves in the direction of the field. 
For large S the shapes develop facets |22j , similar to what 
has been observed for void electromigration [5| • Increas- 
ing the island radius the direction of migration starts to 
deviate from the field direction, and we enter the regime 
of oblique stationary (os) motion (Fig[2J). Since the field 
is aligned with the symmetry axis of the anisotropy, the 
appearance of obliquely moving solutions implies that the 
symmetry of the problem is spontaneously broken. In the 
os regime, pairs of symmetry-related stationary solutions 
coexist; which solution is chosen in a given run depends 
on the initial condition. 

Upon further increase of Rq the obliquely moving 
shapes start to oscillate (Fig|3J). Near the onset of oblique 
oscillatory (oo) motion at radius Rq the oscillation pe- 
riod diverges as T ~ \Rq — R%\~ v with v m 2.5. For 
larger radii higher harmonics of the fundamental oscilla- 
tion period appear and the motion becomes increasingly 
irregular (Fig0J. This characterizes the complex oscil- 
latory (co) regime, which is exemplified in FigEl The 
direction of island motion displays random shifts, which 
seem to be triggered by small fluctuations. This behav- 
ior is typical for large islands, and it is distinct from the 
periodic direction changes seen in the zig-zag (zz) regime 
for moderate sizes and small anisotropics (FigO. 

The true long time behavior for large islands (Rq > 7) 
and large anisotropy (S > 5) could not be pinned down 
unambiguously with our current numerical methods. 
Generally speaking, large islands with small anisotropy 
break up, while for large Rq and S facetted shapes un- 
dergoing irregular motion dominate. 

The example shown in Fig0 provides an important 
clue to the origin of the complex shape evolution. After 
an initial transient lasting until t « 300, the island settles 
down into a shape consisting of a straight upper edge and 
a lower edge which has broken up into a facetted hill-and- 
vallcy structure. The direction of island motion coincides 
with the orientation of the upper, straight edge, as shown 
for a smaller island in the upper panel of Fig0 The key 
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FIG. 4: Time series of the island perimeter showing regular 
and irregular oscillations. From bottom to top the parameters 
are S = 2, R = 5; S = 5, R = 5; and S = 3, R = 8. The 
top panel corresponds to the run shown in Fig[^] 
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FIG. 5: Complex oscillatory motion with S — 3, Ro = 8. 
Light lines show the facet orientations predicted from JHJ, 
dark dashed lines illustrate that the hill-valley structure is 
static in the substrate frame. Consecutive snapshots are 
shifted upwards in time. 

observation is that the hill-and-valley structure on the 
facetted edge does not move in the substrate frame. The 
moving island slides past the static facets, causing the 
shape to oscillate. Around t ~ 900 the roles of the upper 
and lower edges are seen to reverse, and the direction of 
motion changes. 

Quite generally, large islands in the co regime can 
be constructed from four selected facet orientations 
$2, #3; #4- Here B\ and 82 are the possible stable ori- 
entations of the upper island edge, which must satisfy 
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— 7r/2 < 61 < < 6*2 < tt/2. In the case considered here 
(a = 0, n even) the corresponding orientations #3 and #4 
for the lower edge are obtained by reflection at the x-axis, 
6*3 = — 7r — 81 and 6* 4 = it — 8 2 - To form a closed shape, 
at least three orientations must be combined, two of 
which arc those two symmetry-related orientations that 
are closest to the horizontal direction (8 = or ir). In 
Fig0we see a transition from a shape with orientations 
{0i, 83, 84} to a shape with orientations {81,82, 83}. 

The stable facet orientations can be computed along 
the lines of (2?J . The condition v n = for a static shape 
implies that the current in Q is set to a constant j* . 
Using the relation n = dd/ds this can be brought into 
the form 

= -tf*M0)-?**W0)] = -y'ifi\ (3) 

which describes the position 8(s) of a particle moving 
in time s subject to a potential V(8) determined by the 
mobility and the electric field strength. As explained 
in the coexistence of two stable facet orientations 
corresponds to a particle trajectory moving between two 
degenerate potential maxima. To determine the selected 
orientations, j* is tuned until two degenerate maxima 
satisfying the above constraints appear. We have checked 
that this procedure correctly accounts for the facet ori- 
entations observed in the time-dependent calculations 
throughout the relevant region of the phase diagram (see 
Fig|SJ). In general, stable facets can be constructed from 
(J3J) only when the anisotropy is sufficiently large (2^ • For 
n = 6, a = the requirement is S > S c « 1.77. No sta- 
ble facets are found when the number of symmetry axes 
is too small (n < 3); this may explain why we do not see 
oscillatory shape evolution for a threefold anisotropy. 

We close with two remarks concerning future research. 
First, we note that the observed island shapes are quite 
smooth, which implies that the number of circular har- 
monics involved in the shape evolution is small. It thus 
seems promising to attempt a description in terms of a 
low-dimensional dynamical system, in the spirit of |l6| . 
to gain a deeper understanding of the various migration 
modes and the bifurcations connecting them. 

Second, we address the experimental conditions un- 
der which the predictions of this paper could be real- 
ized. As an example, we consider islands on Cu(100), 
for which most materialparameters entering the theory 
are available. Following [8j , we estimate that the electro- 
migration force on an edge atom at a current density of 
10 7 Acm~ 2 is about 400 eV/cm. Together with the exper- 
imentally determined stiffness and mobility f° r 
kinked steps at 300 K, this yields a characteristic length 
of Ie ~ 25 nm, and a time scale ts on the order of sec- 
onds. Thus we expect complex shape dynamics to be 
observable for island radii around 100 nm and on time 
scales of a few hundred seconds. As a first step towards 
a more detailed description of specific surfaces, it would 



be important to identify oscillatory shape evolution in ki- 
netic Monte Carlo simulations of island electromigration. 
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ergy dissipation at surfaces and SFB 611 Singular phe- 
nomena and scaling in mathematical models. 
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